The crossover of deconfinement by tunnelling in FL model 
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We have discussed the tunnelling effect in FL model. The tunnelling coefficient is derived in the 
ffeld configuration space by calculating the transition amplitude using the path integral at SPA 
and the dilute instanton gas approximation. By studying the tunnelling effect between the two 
degenerating vacuums at the critical temperature and chemical potential, we find that the system 
could be deconfined by tunnelling, which will change the first order deconfinement phase transition 
to crossover. The T — fj, phase diagram of deconfinement with both first order phase transition and 
crossover including CEP is presented. 
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I. INTRODUCTION 



It is generally believed that at sufficiently high temperatures and densities there is a QCD phase transition from 
normal nuclear matter to QGP [H, 0]) in which there are theoretically two kinds of phase transitions: chiral phase 

■ transition and deconfinement phase transition. The relation between these two phase transitions is still a puzzle 0] . 
0^ However it seems that the two phase transitions happen coincidently Q. For the chiral phase transition at finite 
Q_j temperature, the quark- antiquark condensate {qRqh) serves as a good order parameter. The order of the phase 
^ _ transition is fairly studied From lattice calculations, for Nf — 2 massless quark flavors, it is a second order 

. phase transition at high temperatures in the chiral limit, while for the non chiral limit it is a crossover. It is difficult 
to make the calculations at finite densities for the lattice, however at small densities one can still derive that the 

■ transition in the non chiral limit is the crossover. At high densities, the chiral phase transition has been studied by 
J> ■ many effective models It is generally regarded that at high densities it is a first order phase transition. Thus 

T-H , in the T — fi phase plane from the first order phase transition at high densities to the crossover at high temperatures 

■ there exists a critical end point(CEP). For deconfinement phase transition, it has not a good order parameter except 
\l [ for infinite quark mass limit, at which the Polyakov loop severs as an ord er p arameter 12, 13]. In the earlier studies 

■ of deconfinement, the bag models, such as MIT model SLAG model [l^ and FriediDerg-Lee(FL) model [ll], had 
^-H ' been often used to investigate the confinement mechanics and the thermodynamics of deconfinement phase transition. 

. With the temperatures and densities increased, the bag will be broken, which results as a first order phase transition 
] ■ of deconfinement. In recent studies the Polyakov loop has been combined into the chiral models, such as Nambu-Jona- 
l ', Lasinio model [13, [3 ^^'^ linear sigma model [19H21I] , which allows to investigate the deconfinement phase transition 
I> ■ within the chiral models. Though the Polyakov loop is not a good order parameter, it still serves as an indicator of a 
transition towards deconfinement. As it is the first order phase transition of deconfinement, the transition will have 
^ latent heat. However, from the recent results of the lattice at high temperatures the deconfinement is also a crossover 
d ' instead of a first order phase transition at the real physical quark mass 0| . 

As the quarks are confined in a nucleon, there is huge binding energy. From hadronic phase to quark phase, 
there should be enormous latent heat, which will result as a first order phase transition. If it is a crossover, the 
transformation between the hadronic phase and the quark phase will be smooth. Thus it is difficult to understand 
where the huge latent heat goes. From the theoretical point of view, what is the physical mechanism lying behind? In 
this paper we wish to use the FL bag model to make a tentative study of the crossover mechanism of the deconfinement. 

The FL model has been widely discussed in past decays [23l - [25| . It has been very successful in describing phe- 
nomenologically the static properties of hadrons and their behaviors at low energy. The model consists of quark 
fields interacting with a phenomenological scalar field a. The a field is introduced to describe the complicated non- 
perturbative features of QCD vacuum. There are two vacuums in this model: one is the perturbative vacuum and 
the other is the non perturbative vacuum or physical vacuum. In the real world the physical vacuum is stable. The 
energy difference between the two vacuums defines the bag constant B. The quarks are confined in a soliton bag. 
It naturally gives a color confinement mechanism in QCD theory. The model has been also extended to finite tem- 
peratures and densities to study deconfinement phase transition [26-30]. In the previous studies, the deconfinement 
phase transition has been often discussed by analyzing the changing of the two vacuums. When the temperature or 
density increases, the energy difference between the two vacuums decreases. Until the two vacuums degenerate, the 
bag constant becomes to zero, and the deconfinement phase transition takes place, which results as a first order phase 
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transition. 

However the tunnelling effect between the two degenerate vacuums are not considered in these studies. It has 
been indicated by T. D. Lee that it is an interesting problem to consider the tunnelling effect between the two 
degenerate vacuums in field configuration space which can eliminate the degenerate vacuums in a finite size microscopic 
system [3l| . Here we think that the tunnelling effect between the two vacuums may be important in deconfinement, 
which could change the first order deconfinement phase transition to crossover. And the tunnelling between the 
degenerate vacuums could be possible physical mechanism for the crossover of the deconfinement. 

The organization of this paper is as follows: in section 2 the FL model is introduced. The field equation of the sigma 
field is treated in homogeneous case. An effective lagrangian is obtained from which a transition amplitude is given 
in a form of the path integral. In section 3, the transition amplitude is evaluated in Euclidean space in stationary 
phase approximation (SPA) and dilute instanton gas approximation. The tunnelling amplitude is further derived. In 
section 4, the tunnelling probability is evaluated for the different effective potentials at different critical temperatures 
and chemical potentials in finite size system in field configuration space. The tunnelling effect is studied in the 
deconfinement phase transition. We indicate that the first order deconfinement phase transition may be replaced by 
a crossover through tunnelling effect. A T — ^ phase diagram of deconfinement with both first order phase transition 
and crossover including CEP is presented. The last section is the summary. 



II. THE EFFECTIVE LAGRANGIAN AND THE TRANSITION AMPLITUDE OF FL MODEL 



where 



We start from the Lagrangian of the FL model, 

£ = ^{ij^d^^ - ga)^p + - C/(a), (1) 

Uia) = ^aa' + ^ba' + ^ca'' + B. (2) 

Tp represents the quark field, and a denotes the phenomenological scalar field, a, b, c, g and B are the constants which 
are generally fitted in with producing the properties of hadrons appropriately at zero temperature. We shift the a 
field as fj (T + cr' where the new a and a' are the vacuum expectation value and the fluctuation of the cr field 
respectively. Then the lagrangian becomes 

£' = i^{i^^d^ - m,)^ + - ^m^'^ - C/(a), (3) 

where 

U{a) = ^aa^ + ^ba^ + ^^ca^ + B. (4) 

lUq — ga and = a + 6(7 + ^ccr^ are the effective masses of the quark and a fields respectively. The interactions 
associated with the fluctuation a' , such as a'^ , a'^ and ipa'Tp, are neglected in the usual tree level approximation. At 
the tree level, one can obtain the fleld equations, 

[zj^d^ - ga)^ = 0, (5) 
d^.d^'a' + mla' ^ 0, (6) 

where (cr''^) and (■(/'V') denote the contributions of the thermal excitations of a and quark flelds respectively. The 
second equality of equation ([T]) means that we define a thermal effective potential. These thermal functions will 
be determined in the later section. Here our discussion will focus on the equation ([7]). cr is the classic field which 
is time and space dependent. There are two physical treatments of a. If one treats the a time independent and 
space dependent which means the static case, the equation ([7]) will be the usual soliton equation which is solved for 
soliton solutions. This case is mostly studied for the confinement mechanics and hadron properties in the previous 
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literature. If one treats the a time dependent and space independent which means the homogeneous case. In this case 
the solution of the equation will be the instanton which is related to the tunnelling effect. Here we are interested in 
tunnelling effect, so we take the homogeneous case and the equation becomes 



da ■ 



(8) 



We can regard the a field as a "space coordinate" of a particle moving in one dimension in a potential 14// (cr) . By 
this analog we can further obtain the effective Lagrangian density which describes the motion of the particle, 



1 f da 



^^ff 2 V dt 
For the homogeneous case the effective Lagrangian is 



Vesf{a). 



L 



d'xC 



1 fda 

2 I 



Veff{a) 



(9) 



(10) 



in which f2 is the space volume. The effective Lagrangian could be quantized by the canonical quantization. The 
corresponding quantized Hamiltonian in quantum theory is 



^=|^ + f^K//(a), 



(11) 



where p = —ihd/da is the canonical momentum operator. Consequently the quantum transition amplitude for a 
particle to start at "position" a — aa &t time t = 0, and end up at cr = (Tft at t = T is given by 



K{ab,T;aa,Q) =< ab\exp{-iHT/h)\aa > 
By the path integral the transition amplitude can be also written as 



where 



S[a{t)]^ f dtLeff= [ dm 
Jo Jo 

is the action and / D[a{t)] is an integral over all path. 



Kia,,T;aa,0)= / D[ait)]exp{i/hS[ait)]}, 



1 fdaV 

2U) -^^^^('^^ 



(12) 



(13) 



(14) 



III. INSTANTON AND TUNNELLING EFFECT 



In order to evaluate the transition amplitude, it is more convenient to do the calculation in a Euclidean space. The 
corresponding Euclidean transition amplitude denotes as 



KEiab,T/2;aa,-T/2) = / DE[aiT')]exp[~l/hSEiaiT')], 



where 



SE[air')]^ f'\r'^ \\ (|,)%K//(.) 



is the Euclidean action, t = iT and r' = it. One can also obtain the Euclidean equation of motion, 

d^a ^ dV^ffja) 
dr'^ 



da 



(15) 



(16) 



(17) 



In the following discussion we take r — > 00 and we suppose V^f / (a) has the configuration schematically sketched in 
FiglH The potential has two minima which are degenerate. The equation (|17p could be solved by the boundary 



cr 



FIG. 1: The schematically sketched thermal effective potential as a function of cr. 
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FIG. 2: A single instanton solution in the FL model. 



condition 

t' — —CO, a = Ua', t' — oo, a — Gb- (18) 

The solution is an instanton which is schematically sketched in FiglJJ This instanton solution aci{T') is a classical 
path which is an extremum of the action SeIc^t')]. One could expand the action SeIc^t')] about the classical path 
a- = (Tci{t'), 

SeHt')] = SEW,i{r')] + - / dT'y{T')OE{r'}y{r') + 0{y'), (19) 



where 



and 
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y(r') = a(r')-a,Kr'), (20) 

In the spirit of the stationary phase approximation (SPA), the major contribution to the transition amplitude (jisp 
will come from the vicinity of the classical path cr(r') = adir'). In that case, to leading approximation one can 
neglect the cubic and higher terms in ?/(t') in the Taylor series The SPA as applied to the path integral in the 

transition amplitude gives 

KEi<7b, r/2; -r/2) = e-^°/''SB(r) [DetOEir')]-^ , (22) 
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where 5*0 = S'^ftTciC'r')] and Be{t) is the measure factor in functional integral. Det is the determinant of the operator 
Oe(t') which is most conveniently obtained from its eigenvalues, but the operator will have a zero eigenvalue because 
the action is translation invariant in the r' variable. Formally this will cause a divergence in (|22|) . However this zero 
mode could be separated from the determinant. For a detail discussion of this treatment one could refer to [s^]. After 
separating the zero mode the result is 

KE{ab,T/2-aa,-T/2) = e-^^''' JtBe{t) [Det' Oe{t')]~^ , (23) 

where J = \/Sq and Det' denotes the determinant without the zero mode. In order to calculate this determinant, 
one could be referred to the case of harmonic oscillator. If the particle oscillates in the vicinity of the "position" aa, 
it could be regarded as a harmonic oscillator. So is it in the vicinity of the "position" ab- We suppose the angular 
frequency of the harmonic oscillator is w. The determinant can be written as 

k, (24) 

where fc is a constant independent of r as r oo. For the harmonic oscillator, the determinant could be evaluated 
at T — ^ oo and the result is 

BE{r) 

Hence the amplitude (l23l) becomes 

KE{,ab,T/2-aa,-T/2)^e-^°''\JkT(^^~^e'^. (26) 

From the above derivation we obtain the result of the transition amplitude by considering single instanton. From 
the notion of tunnelling, this means the particle traverses the potential barrier once. However the particle could also 
traverse the potential barrier several times, that means there are also multi-instanton solution which configuration 
is shown in Figl3]that will have the correct boundary condition for our problem, namely a{—oo) = tJa, a{oo) = ab- 
The instantons should be widely separated and is so called the instanton gas. For our interest of tunnelling effect. 
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(25) 



FIG. 3: Multi-instanton solution in the FL model. 

the contributions from all these multi-instanton solutions to the transition amplitude should be summed up. The 
summation could be evaluated under the dilute instanton-gas approximation [32| and the result is 

1 



KE{<Jb, r/2; (7^, -r/2) = i^—j ' e~^sinh i^Jkre-'^^/'^j . (27) 
By equation (|12|) the transition amplitude can be also written as 

< ab\exp{-HT/h)\aa >= ' e-'^sinh (jfere"^"/'*) . (28) 
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By the derivation in a similar way to the above derivation one can also obtain three kinds of transition amplitudes in 
the following, 



C (Ta\exp{-HT/h)\aa >= (1^) ' e^'^cosh (^Jkre'^"^'^^ , 
< ab\exp{-HT/h)\ab >= ' e-'^cosh {^Jkre'^"'^ 



If we define the states 



k± >= ^(ka >±|ct6 >) 



then from equations ([251 — f?T|) we have 



<a±\exp{-HT/h)\a± >= ( ^) %-^±^*^— o/«^ 



From the above expression one could read out the energy 

1 . 



(29) 
(30) 
(31) 

(32) 

(33) 
(34) 



from which one could see that the states |cr± > are the eigen-states of H with the eigenvalues E±. The first term 
in equation p4p is the energy without tunnelling while the second term is the energy correction coming from the 
contribution of the tunnelling through the potential barrier. The factor e~^°^^ reveals the exponential dependence 
on the "barrier strength" . Notice that 



/oo 
'OO 



(35) 



in which CTci(T') is the classical path and will obey 



1 f da 



O =^e//(a). 



(36) 



2 \dT' 

The above equation can be derived from the equation p7)) by integration. Substituting equation ((36)) into equation ([35 
and changing the integration variable we obtain 



5*0 = /" daflJ2Veffia)da. 



(37) 



This quantity represents the strength of the potential barrier between ctq and Cf,. It is clear that the tunnelling 
amplitude is 



M — exp 



Sq 



exp 



f ' daJ2Veffi<j)da 



(38) 



The result is in accord with conventional WKB treatments of quantum mechanics in the position configuration space, 
but one should notice that the result here is in the field configuration space. 



IV. CROSSOVER OF DECONFINEMENT BY TUNNELLING 



By the tunnelling amplitude one can further obtain the tunnelling probability 

2VL 



ry = |M| = exp 



2Veff{a)da 



(39) 
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It is clear that the tunnelhng probabihty is dependent on the potential Veff{(j) and the volume ft. The effective 
potential V^ff should be determined in finite volume J7. However first we will determine the effective potential in 
the infinite volume limit. Recall equation ([7|), the thermal excitations of a and quark fields could be evaluated by 
standard method of finite temperature field theory and the results are 



in which /3 is the inverse temperature, fi is the chemical potential, and 7 is a degenerate factor, 7 = 2{spin) x 

2{flavor) x 3(coZor). E^r = ■\/p-^ + and Eq = y^p^ + rn^. For the values of the model parameters a, &, c and 5, 

there are different choices. We have taken one set of values as a = 17.7/to^^, h = — 1457.4/to^^, c ~ 20000, g = 12.16, 
which have been often used in the literature 28]. The effective mass of a field is fixed at to^ — SSOA/eV. One should 
notice that the second equality of equation ([7]) means that we define a thermal effective potential. By equation (|40|) 
and (|4ip the thermal effective potential in the infinite volume limit is 

= U{a) + ^ / ln(l - e-^^^) " ^ / [ln(l + er'^^^^-^^) + ln(l + e-^(^'+-))] . (42) 

Next we will consider the finite volume effect. In the case of a finite system of linear size L the momentum integral 
in the effective potential will be replaced by a sum, which means 

Veffia) = Uia) + ^ E Ml - e-^^-) ~ j^Y. [^1 + e-'^(^-^)) + ln(l + e-^(^'+''))] . (43) 

For a cubic system with the volume Vt = one should consider the boundary conditions. In the paper [s^], it is 
indicated that one is free to choose the periodic boundary condition (PBC) or the anti-periodic boundary condition 
(APC), which will lead to different discretization in the summation. In the limit of L — > 00, the two different 
discretizations all go to the same continuum limit. However in the realistic evaluation we find that APC will be more 
consistent with the intrinsic statistic property of the fermion, while PBC is more proper for the boson. Thus we choose 
the PBC on the summation relating to the sigma field, while the APC on the summation relating to the quark field. 
That means in equation (143^ for the summation in the second term the discretized momentum is ki = 2iTmi/L^ while 
for the summation in the third term the discretized momentum is U — 7r(2rti + l-)/L, where rrn and n.i are integers 
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nii, fii = ■ ■ ■ , —2, — 1, 0, 1, 2, • • • and i denotes for x, ?/, z. For E^ and Ei we have Ek — \J'^ + and Ei ~ + m^, 

in which = fc^ + fc^ + k\ and P = + + 'z- To determine the system size L, we mention that the charge radius 
of a nucleon is about Tc ~ 0.8/to. For a nucleon in a cubic box with linear size L, one can take the value of the size 
L « 2rc ~ 1.6 fm (3^ . Then the summations in equation (j43p could be numerical evaluated at the given temperature 
and chemical potential. 

Now we are in a position to discuss the effective potential at different temperatures and chemical potentials in a 
finite size system. At zero temperature and chemical potential, the Veffic) is just the U{a), which has two minima 
at (Ta = and ct;, 7^ corresponding to the two vacuums. The CTb 7^ is the absolute minimum which is the physical 
vacuum. When the temperature or chemical potential increases, by evaluating the effective potential from equation 
(US]) we find that the energy difference between the two vacuums decreases. Until the two vacuums degenerate that 
the bag constant B becomes to zero, the deconfinement phase transition takes place, which determines the critical 
temperature and the critical chemical potential. In Fig|4]we show the thermal effective potentials at different critical 
temperatures and critical chemical potentials in the finite size L. It is clear that the phase transition here is first order. 
However one should notice the tunnelling effect between the two degenerating vacuums. When the tunnelling takes 
place, the physical vacuum will be smoothly transformed to the perturbative vacuum without any energy changing, 
and there will be no latent heat as a result. The first order phase transition will be replaced by a smooth crossover. 
The tunnelling probability depends on the effective potential. From FigUlwe can see the potential barriers become 
fiatter and narrower from the top one to the bottom one, which will result a smaller tunnelling probability by a 
qualitative analysis. In order to make it more accurate we will evaluate the tunnelling probability. 

In the calculation of tunnelling probability, we take S = 1 and the system volume \s ^ — ^ 4.1 fm^. From 
equation we could evaluate the thermal effective potential at the critical temperature and chemical potential. 
By equation p9p the tunnelling probability could be further determined numerically. The tunnelling probabilities 
for different critical temperatures and chemical potentials are shown in Table ID From the table we could find that 
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-0.05 0.00 0.05 0.10 0.15 0.20 
o-(fm"') 



FIG. 4: The thermal elfective potentials at the dilferent critical temperatures and chemical potentials. From top to bottom the 
critical temperatures and chemical potentials are: Td = 5MeV, fid = 668MeV; Tc2 = 60MeV, fic2 = 480MeV; Tcs = 12QMeV 
Hc3 = 239MeV; Tc4 = 153MeV fid = OMeV; 



TABLE I: The tunnelling probabilities at the different critical temperatures and chemical potentials. 



Te(MeV) iJ.c{MeV) rj_ 

5 668 49.5% 

20 628 56.0% 

40 557 63.8% 

60 480 68.1% 

90 362 72.6% 

120 239 75.1% 

153 77.1% 



the tunnelling probability increases with the critical temperature increasing (or with the critical chemical potential 
decreasing). That is to say at the high critical temperatures and the low critical chemical potentials along the phase 
boundary of deconfinement it is more probable for the tunnelling to occur than that at the low critical temperatures 
and the high critical chemical potentials. From Table U we can set a threshold value of tunnelling probability as 

r]c w 75.1%. (44) 

When r] > rjc, we think that the tunnelling effect will dominate, and the transition from hadronic matter to quark 
matter is a crossover by tunnelling. When rj < rjc, the tunnelling effect will be weaker, and the transition will be 
dominated by the first order phase transition, rj^ is a critical point which value may not be the one we choose 
here, however there should exist a certain value of the critical point that separates the tunnelling process and the 
first order phase transition. This critical point is regarded as a critical end point of the first order phase transition. 
Therefore we can plot the T — fi phase diagram of deconfinement phase transition as shown in FiglS) It is clear that 
the deconfinement phase transition is first order at the low critical temperatures and high critical chemical potentials 
while crossover at the high critical temperatures and low critical chemical potentials. From the first order transition 
to the crossover there exists a CEP. 



V. SUMMARY 



In this paper we have discussed the crossover of deconfinement by tunnelling in FL model. It is indicated that the 
system could be deconfined by tunnelling process which will lead to a crossover instead of a first order deconfinement 
phase transition. By the method of path integral we have derived the tunnelling amplitude in the field configuration 
space and obtain the tunnelling probability which is accord with the conventional WKB treatment of quantum 
mechanics. The tunnelling probabilities for the different thermal effective potentials in the finite volume in field 
configuration space at the different critical temperatures and chemical potentials are evaluated. We find that the 
tunnelling is more likely to take place at the high critical temperatures and low critical chemical potentials. A 
threshold value of tunnelling probability is set. Above that value the system is deconfined by tunnelling, in other 
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FIG. 5: T — fx phase diagram of deconfinement with CEP in the FL model. 



words it is a crossover of deconfinement. Below that value it is a first order deconfinement phase transition. The 

T — II phase diagram of deconfinement with both crossover and first order phase transition including the CEP is given. 
However it should be reminded that in this paper we only give one possible physical mechanism which could lead to 
the crossover of deconfinement, and the corresponding numerical results on the phase diagram are model dependent. 
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